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SUMMARY 


Flight tests with the NASA Terminal Configured Vehicle (TCV) Boeing 737 
were made to investigate the use of back-azimuth signals to estimate airplane 
position. Position estimates were established by three update methods; a 
combination of microwave landing system (MLS) back-azimuth and distance- 
measuring-equipment (EME) signals, MLS back-azimuth signals alone, and dual DME 
signals. Position-estimate error was obtained by comparison of the position 
estimate from the airplane navigation system with the tracking information from 
a ground-based radar system. 

The results of these tests showed that the most accurate position esti- 
mates were obtained with a combination of back-azimuth and DME signals. The 
next most accurate estimates were made with only back-azimuth signals. Of the 
three update modes tested, the least accurate estimates were made with dual 
EME signals. 

Analysis of the position-estimate error data showed that the component of 
error due to the EME signal input was larger than the component due to the MLS 
back-azimuth signal input. The use of a precision DME signal input would be 
expected to reduce position-estimate error. 


INTRODUCTION 

The NASA Terminal Configured Vehicle (TCV) previously demonstrated the 
ability to navigate along a curved path approach to landing (ref. 1) with an 
onboard navigation system utilizing position information from the time- 
referenced scanning-beam (TRSB) microwave landing system (MLS) . During 
operations in the MLS environment, azimuth and elevation-angle information was 
obtained from two scanning beams, and distance information was provided by MLS 
precision distance-measuring equijment (PDME) for use in the onboard navigation 
and automatic landing systems. The TRSB MLS provides azimuth information 
within 60° of each side of the runway center line, elevation-angle information 
up to 20°, and FEME range information up to 37 040 m (20 n. mi.) from the MLS 
antennas . 

Missed-approach navigation capability was available from the MLS system 
in the form of back-azimuth and PDME information. Hence, to fully utilize the 
TRSB MLS, algorithms using back -azimuth information were incorporated into the 
TCV airplane navigation computer for use in generating missed-approach guidance. 
Flight tests then were conducted to determine the navigation accuracy from MLS 
back-azimuth signals. 

The purpose of this report is to describe the equations and logic used to 
generate a navigation position estimate in the MLS back-azimuth signal environ- 
ment and to document the error in the navigation position estimate. The equa- 
tions described are used to calculate position difference components from 
(1 ) MLS back-azimuth bearing and range information from arbitrarily located 



DME stations, and (2) MLS back-azimuth bearing information alone. A summary 
of the TCV position-estimate update process is also described. The navigation 
position-estimate error calculated from flight data and radar tracking infor- 
mation is analyzed to determine the relative errors due to DME and back-azimuth 
signal inputs. The position-estimate error data using the MLS inputs are also 
compared with error data obtained during dual DME updates, the primary position- 
estimate mode. 
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SYMBOLS 

distance between the DME and MLS back-azimuth antenna, n. mi. 
east component of A, n. mi. 
north component of A, n. mi. 

distance from DME to airplane corrected for slant range, n. mi. 
slant range distance from DME to airplane, n. mi. 
position difference vector 

component of position difference vector DP perpendicular to the 
runway center line 

magnitude of vector DPp, n. mi. 

magnitude of the component of position-estimate error perpendicular 
to the back-azimuth bearing on which the airplane is located, m 

magnitude of the component of position-estimate error parallel to 
the back-azimuth bearing on which the airplane is located, m 

azimuth bearing-error angle, deg 

ellipticity constant, 0.003367 

altitude of airplane above mean sea level, m 

altitude of DME antenna above mean sea level, m 

altitude of MLS back-azimuth antenna above mean sea level, m 

unit coordinate vectors 

position update gain 

velocity update gain, sec“^ 

outer radial limit of back-azimuth volumetric check, m (n. mi.) 
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inner radial limit of back-azimuth volumetric check, m 

axes of orthogonal coordinate system oriented toward true north 

angle formed by the vector and a line between the airplane and 

DME antenna, deg 

range from the back-azimuth antenna to the airplane, m 

radius of Earth, m 

meridional radius of curvature, m 

normal radius of curvature, m 

unit vector perpendicular to the runway center line 

north, east components of navigation system velocity estimate, knots 

north, east components of inertial ground speed, knots 

axes of orthogonal coordinate system oriented along the runway 
center line 

->• 

coordinates of vector Zq transformed into the X',Y' coordinate 
system 

vector of airplane estimated position from the MLS back-azimuth 
antenna 

magnitude of Zg, n. mi. 

-V 

north, east components of Zg, n. mi. 

vector of airplane measured position from the MLS back-azimuth 
antenna 

magnitude of Zj^, n. mi. 

north, east components of n. mi. 

vector of airplane position estimated radially along the measured 
azimuth angle 

magnitude of position vector Z^, n. mi. 

north, east components of position vector Zj., n. mi. 

angle formed at the back-azimuth antenna by the DME antenna and 
the measured airplane position, deg 


APjj, APe 


north, east components of position-estimate error, n. mi. 



At 


change in time 


AVjgjAVg north, east components of system velocity update, knots 

A(j),AX latitude, longitude update estimates 

f| back-azimuth angle relative to the runway center line, deg 

)j relative angle between the DME and the MLS back-azimuth antenna, deg 

'I’DME'^DME latitude, longitude of DME antenna location, deg 

<})g,Xe latitude, longitude of airplane position estimate, deg 

4)q,Xq latitude, longitude of back-azimuth antenna location, deg 

\l^j. runway heading to true north, deg 

vertical angular limit of back-azimuth volumetric check, deg 
(jj lateral angular limit of back-azimuth volumetric check, deg 

Subscripts: 
t at time t 

t-1 at iteration time previous to time t 

Abbreviations: 

DME distance-measuring equipment 

IBD navigation update mode: inertial velocity, back azimuth and DME 

IBX navigation update mode: inertial velocity, back azimuth alone 

IDD navigation update mode: inertial velocity, dual DME 

MLS microwave landing system 

NCU navigation computer unit 

PDME precision distance-measuring equipment 

VOR very-high-frequency omnidirectional-range ratio 
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AIRPLANE, EXPERIMENTAL SYSTEMS, AND EQUIPMENT 

Airplane 


The NASA test airplane is the Boeing 737-1 00 twin-jet transport airplane 
shown in figure 1 . The TCV airplane is used as a research vehicle with 
separate experimental navigation, guidance, flight control, and display 
systems located in a separate research cockpit. All of the normal flight 
systems (flight control, navigation, pressurization, etc.) have been retained 
in the conventional cockpit in a normal, functional state. This arrangement 
allows changes to be made to any of the experimental systems while retaining 
the standard operational features of the airplane. 


Experimental Systems 

The experimental research flight systems consist of a triplex digital 
flight control computer system, an electronic cathode-ray-tube (CRT) display 
system, and a digital navigation and guidance system (ref. 2) integrated into 
a separate, two-man-crew, research cockpit. Figure 2 is a simplified func- 
tional block diagram of the navigation, guidance, and control process during 
automatic path tracking (ref. 3) . The navigation and guidance calculations 
are performed in a single digital navigation computer. Various navigation 
sensor signals (including inertial-navigation-system velocities and accel- 
erations, true airspeed, magnetic heading, and MLS back-azimuth, VOR, and 
DME radio signals) are used in the navigation computer. This computer esti- 
mates the airplane position based on combinations of the sensor inputs. 
Horizontal path, vertical path, and thrust commands based on the estimated 
position, velocity, and path tracking errors of the airplane are computed 
and transferred to the flight control computer system 20 times per second. 

The flight control system then commands the flight control-surface servos. 


Airborne Data Acquisition System 

A wide-band magnetic tape recorder onboard the airplane recorded data at 
a rate of 40 samples per second. These data included 93 parameters describing 
the airplane configuration, attitude, and control-surface activity. Thirty-two 
additional channels were used for recording parameters, including the latitude 
and longitude of the navigation position estimate which was calculated by the 
navigation computer. Video recordings of the electronic attitude-director and 
electronic horizontal-situation displays were also made throughout the flights. 


Radar Tracking Facilities 

Radar tracking of the airplane was provided by the Federal Aviation 
Administration Extended Area Instrumentation Radar (EAIR) Facility at the 
National Aviation Facilities Experimental Center (NAFEC) in Atlantic City, 
New Jersey. The tracking radar is a precision, C-band, instrumentation 
radar system which was operated in a secondary (beacon-tracking) mode during 
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these flight tests. Slant range, azimuth angle, and elevation angle data 
were recorded at a 10-Hz sample rate on magnetic tape. All airborne and 
ground-based recorded data were time correlated for post-flight processing 
and analysis. Post-flight processing of the range, azimuth, and elevation 
data consisted of conversion to latitude, longitude, and altitude. 

Radar accuracy, root mean square (rms) , is stated (ref. 4) as 0.15 mrad, 
5.6 m at 37 040 m (20 n. mi.) in azimuth and elevation and 6 m in range. 


MLS and EME 

The MLS used during these flight tests was a test bed system oriented 
for approaches to runway 31 at NAFEC. This system was a time-referenced 
scanning-beam system with a front-azimuth beam width of 1°, a lateral sweep 
angle of +60°, and nominal vertical coverage up to 20°. The glide slope 
beam was 1° wide and had an azimuthal coverage of approximately ±60°. An 
MLS PDME with omnidirectional coverage of approximately 37 040 m (20 n. mi.) 
was used in conjunction with the front-azimuth signals. The back-azimuth 
beam was 3° wide and had a lateral sweep angle of ±40°. Vertical coverage 
was between 1° and 20°. 

The Atlantic City DME, normally used for en route navigation, was used 
during the back-azimuth flight testing. This EME is located at the airport 
1671 m past the back-azimuth antenna and 229 m left of the runway center 
line (fig . 3) . 

When the navigation system updated its position estimate with dual DME, 
the Atlantic City DME was always used in conjunction with other DME in the 
area. The navigation system automatically selected, tuned, and checked the 
validity of the other EME on the basis of strength of signal and station 
location geometric characteristics so that good position estimates resulted. 
These geometric characteristics are described in reference 2. 


Switching Logic of Airborne MLS Receiver Antenna 

The TCV airplane was equipped with two MLS receiver antennas to ensure 
that the greatest possible signal strength was delivered to the MLS receiver 
regardless of the airplane direction of flight, its attitude, or its position 
within the MLS signal coverage. The forward antenna was located on top of the 
fuselage just above the forward cockpit. The rear antenna was mounted on 
the bottom of the fuselage below the airplane tail. The MLS airborne receiver 
contained the logic to determine which antenna was providing the greatest 
signal strength. The receiver sent a pair of discrete logic signals to an 
external antenna switch to select the appropriate antenna. 

Additional hardware logic in the flight control computers forced selection 
of the forward antenna while operating in the MLS front-course signal coverage, 
regardless of the relative signal strengths at the forward and aft antennas. 
This additional antenna-selection logic was required because the flight control 
computers did not have the necessary software for antenna switching. 
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EXPERIMENT DESIGN AND PROCEDURES 


The objective of this flight experiment vras to determine the accuracy 
of the navigation position estimate when updated with dual DME signals, with 
MLS back-azimuth signals, and with a combination of MLS back-azimuth and 
IME signals. The flight and radar data obtained during these tests were used 
to compare the accuracies of these update modes and to determine the relative 
errors due to DME and back-azimuth signal inputs. 

The position-estimate error was determined by calculating the vector 
between the estimated position of the airplane and its position as tracked 
by ground-based radar (ref. 5) . The magnitude of this error vector is the 
position-estimate error . The mean and standard deviation of the navigation 
position-estimate error was calculated and used for comparison of the three 
update modes . 

The NASA TCV airplane was flown along three, three-dimensional flight 
paths (fig. 3), each stored in the navigation data base in the navigation 
computer. Path 1 was a 5556-m (3-n. mi.) downwind leg with a 180°, descending 
left turn to a 5556-m (3-n. mi.) final approach to the runway. This path was 
in the MLS front-course signal coverage and allowed the navigation computer 
to update its position estimate with the front-course signals just before 
each test run. 

Paths 2 and 3 provided experience with guidance in the area of back- 
azimuth signal coverage. These paths overlaid the last 5556 m (3 n. mi.) 
of path 1 so that lateral guidance would be continuous during the transition 
from front-course to back-course guidance. Path 2 extended 18 520 m (10 n. mi.) 
from the runway threshold, straight out along the runway center line. Path 3 
followed the runway center line for approximately 6482 m (3.5 n. mi.), then 
made a 45° left turn, and continued an additional 21 298 m (11.5 n. mi.). 

Each test run was started with the pilots flying an instrument approach 
to the runway via path 1 . Instead of landing, a missed approach was executed, 
and the airplane was flown in the MLS back-azimuth coverage along either path 2 
or path 3. Upon reaching the end of either of these paths, the airplane was 
returned to the beginning of path 1 for another test run. 

The radar was calibrated with the airplane on a known ground test point. 
The radar tracked the airplane continuously during the flight from take-off 
to landing. 


DESCRIPTION OF NAVIGATION POSITION ESTIMATE 
Position Difference 

The TCV navigation computer was programmed to select and tune two appro- 
priate DME and/or VOR stations in the vicinity of the airplane. The distance 
and/or azimuth information received from these stations was used to update 
the previous position estimate. This position difference was divided into 
north (APjg^t) and east (APg^^) components to be used in determining new posi- 


tion and velocity estimates. The dual EME update mode, the primary position- 
estimate difference mode, was displayed as IDD on the electronic map displays 
for each of the flight crews. 

Other sources of navigation information could also be used to determine 
position differences. While operating in the MLS back-azimuth signal cover- 
age, these position-update components were calculated with back-azimuth and 
DME information. This update mode was shown on the pilots' electronic map 
display as IBD (inertial velocity, back azimuth, DME) . If no DME was available, 
the estimate mode was IBX (inertial velocity, back azimuth alone) . In subsequent 
sections, the equations used in the IBD and IBX update modes are presented. 

The equations used in the IDD update mode are presented in reference 2. 

In the TCV navigation system the same position-estimate calculations, 
described in the next section, were used regardless of the position-update 
mode used to determine the north and east components. 


Position Estimate 

The first step in the position-estimate process was to develop a velocity- 
estimate update from the position update: 


^^E,t = + ^2 ^^E,t 

An updated velocity estimate was then made by summing the velocity update 
with ground speed obtained from the inertial navigation system: 


%,t - ^%,t + ^N,t 
^E,t = ^^E,t ^E,t 


A position update in terms of latitude and 
the system velocity and the position-difference 


longitude was obtained from 
components as follows: 


A4>t 


VN,t At + Kt APN,t 

*^M,t 


AXt 


^E,t At + Ki APg^t 

*^N,t 
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This latitude and longitude position update was based on an oblate-spheroid 
Earth model (ref. 2) by using the following radii of curvature in the meridional 
(north/south) and normal (east/west) directions: 


*^M,t = ha/p,t + ^£(1 - 2F + 3F sin2 (j)e,t-T) 

^^N,t = [ha/p,t + rsO + F sin2 cos 4)e,t-l 


The updated position estimate was found by summing the previous position 
estimate with the position- update terms; 


‘5’e,t = 'l’e,t-1 + ^'t’t 

^e,t = ^e,t-l + ^^t 

Dual DME (IDD) Position Estimate 

Dual DME updating is the normal position update mode. The equations used 
to determine position difference while in this mode are given in reference 2, 


Back-Azimuth/DME (IBD) Position Estimate 

General solution .- Figure 4 shows the geometry of the MLS back-azimuth 
antenna, the DME antenna, the position estimate of the airplane, its position 
measured with back-azimuth and DME information, and the runway with an extended 
center line. An orthogonal coordinate system, with its origin placed on the 
back-azimuth antenna is oriented with respect to true north. Since the DME 
could be selected either manually or automatically, the position-update equa- 
tions were derived so that the DME was not required to be colocated with the 
back-azimuth antenna. 

By determining the relative geometry between the back-azimuth antenna, the 
DME antenna, and the position of the airplane measured with the back-azimuth and 

DME information, the vector between the origin and the measured position 

of the airplane was found. The vector Zg ^-1 between the origin and the pre- 
viously estimated position of the airplane could be found directly from their 
known latitudes and longitudes. The position difference was then found by sub- 

tracting from . 
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Calculation of the airplane position vector Figure 4 shows that 

Zm, t side of the triangle formed by the origin (back-azimuth antenna), 

the measured position of the airplane, and the DME antenna. Known quantities 

used to determine Zjj, ^ in this triangle include the runway heading the 

slant range of the DME D^, the back-azimuth angle relative to the runway center 
line ri, the location (latitude and longitude) of the DME antenna ^EME' 

and the back-azimuth antenna location <1 )q,Xq. 

Figure 5 shows the angular geometry and distance A between the origin 
and the DME antenna. The length A and the angle between the DME antenna 
from the north axis y remained constant for a particular DME and were cal- 
culated only once. If a different DME was tuned, then A and y were 
recalculated. The length of A in nautical miles was determined by vec- 
tor i ally summing its components, a and b; 


I 'I’o 

® ~ (^DME ~ ^o) ^02 I 

b = ('t’oME ~ 'f’o^ ( 60 ) 




The angle y was found by 


y = tan ^ 



(0 S y S 21T) 


The angle a, formed by side A and vector may vary continuously 

and was calculated 20 times per second. It ranged between 0 and tt, inclusive, 
and was found from 

at = l 2 iT - y + (ilJj. - Tit) I 


If this calculated value of at was greater than tt, then 


«t 


2tt - [l27r - y + (Tpj. 



The magnitude of the DME reading Dt, measured in the airplane, was the 
slant range distance between the ground-based DME antenna and the airplane. 


10 


This distance was slant-range-corrected to determine the ground distance 
between the airplane and DUE 


Dt = Di. sin 


cos 


-1 


- h' 


CME) 


Angle P, formed 
and was calculated 20 


by the vector side D^f could vary continuously 

times per second. From the relation. 


A °t 

sin sin Oj. 


angle P|- was determined from 


P^ = sin ^ 



sin a± 


Comparison of the square of side A with the sum of the square of side D^. 
and the square of the magnitude of the vector (estimated airplane 

position vector) determined whether angle P^- was obtuse or acute. Hence, 


sin~^ ( — sin a* 


Pt=< 


TT - Sin 


/ A 

-1 


I — Sin a*. 

\^t ; 


Lf a 2 ^ Dt2 + Zg^t-1 


if a2 > 


The magnitude of the estimated airplane position vector was used as an 

approximation for The magnitude of Zg^^._-| was 


i,t-1 = + (2g^E,t-l) 


n 



in which 


2e,N,t-l = ~ (l)o)(60) 

Ao + 

The actual magnitude of Zjj,^ t was found from 
^m,t “ ^ ^‘t ■*■ ^t ^t 

The north and east components of were calculated from the angle - ri^. 

->■ 

between the north axis and 

Zm,N,t = Zm,t cos - Ht) 

Zm,E,t = ^m,t sin (i^r - Ht) 

Calculation of the position difference in north and east components ^N,t' 
APg^l- — IBD update mode.- The position difference in north and east components 

was now found by subtracting the north component of from the north 

“>■ ->■ 
component of '^xa,t subtracting the east component of from the 

->■ 

east component of 

^N,t = - Zij,^N,t 

^E,t = Zg^E,t-l ” 2m,E,t 

These position difference components were then used directly in the navigation 
position-estimate algorithms. 
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Back -Azimuth Only (IBX) Position Estimate 


General solution .- In the event that a valid DME signal could not be 
obtained, the navigation computer utilized the back-azimuth signal and inertial 
velocity to determine a position estimate. Figure 6 shows the geometry of the 
back-azimuth antenna, the position estimate of the airplane, its estimated 
position on a measured back-azimuth radial, and the runway with an extended 
center line. An orthogonal coordinate system, with its origin at the back- 
azimuth antenna (4 )q,Xq), is oriented in a true-north direction. 

Since no DME information was available, no radio position difference could 
be developed along a radial from the back -azimuth antenna. However, inertial 
velocity was utilized in the radial direction to supply inputs for a new posi- 
tion estimate in the position-estimate algorithms. Radio position difference 
in the IBX update mode was limited to a direction perpendicular to the runway 
center line. 

North and east position estimates were found in the following manner. A 

position-difference vector DP^- was found by subtracting the estimated position 

vector from the estimated position vector ,t measured back- 

azimuth angle. The component of DP^. perpendicular to the runway was found 
and broken into north and east components. These components were used in the 
navigation position-estimate algorithms. This process was repeated 20 times 
per second . 

Calculation of DP^..- To determine DPj- it was necessary to calculate the 

position-estimate vector Zg ^-1 north and east components from the latitudes 

and longitudes of the last position estimate and of the back-azimuth antenna 
location; 


Ze,N,t-1 - (60) (cj)g^t-l ■ 


2e,E,t-1 


(60)(Xe,t-l 


Xq) cos 




A vector Zj. ^ of the estimated position of the airplane along a measured 
back-azimuth radial was determined in the following manner. Since no radio 
updates could be obtained along a radial from the back-azimuth antenna, it was 
assumed that the estimated radial distance from the origin was correct. Hence, 
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the vector lengths of Zj.^^ and were the same (directions could differ 

' > 

to obtain DP|.) . The length of Zj-^f. and was 



The north and east components of were found from the angle vIJj- - 

between the vector Zj-^t north axis: 


Zr,N,t = Zr,t cos - nt) 

^r,E,t “ ^r,t sin (4^^ “ ^t^ 

^r,t ~ (^r,E,t^^ (^r,N,t^i 

DP^ was found vector ially: 

°^t “ ^r,t ■ ^e,t-l 

Calculation of the position difference APjj^^/ '^^E,t “ update mode.- 
The magnitude of the component of DP^ perpendicular to the runway center line 
was obtained by the dot product of DP^. and a unit vector u perpendicular to 
the runway center line- The unit vector u, shown in figure 6, was 

/N /\ -'N 

u = -COS (4’j-)i + sin (4'r)j 
The magnitude of DPp^t thus 

DPp,t = DPt • ^ 

= -(Zr,E,t - Ze,E,t-l) cos ipj. + (Zj-,N,t " Ze,N,t-l) sin <})r 


This provided the magnitude of the position difference. The north and east 
components of the position difference were 
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APN,t = Sin iij. 


COS 'l^r 


These position-difference components were used directly in the navigation 
position-estimate algorithms. 


Back-Azimuth/Airplane Position Validity Check 

Since random DME selection, automatic frequency tuning, and other means 
of automatic software control were utilized in the navigation computer, a check 
had to be made to ensure that the navigation computer was using the appropriate 
navigation data for the area in which the airplane was flying. Obviously, if 
improper navigation data were being utilized, position estimates based on those 
data would be inaccurate and should not be used. 

To preclude the use of the wrong MLS back-azimuth navigation data, a volu- 
metric geometric check was made to determine whether the position estimate of 
the airplane was within the MLS back-azimuth boundaries. If the airplane was 
not within these boundaries, then MLS updating was inhibited. These boundaries 
included lateral and radial limits of coverage as shown in figure 7. The 
lateral limit of the back-azimuth angle of coverage is ±to from the runway 
center line. Radial limits required that the airplane be within a distance L 
in meters (n. mi) of the back-azimuth antenna, but not closer than a distance M 
in meters (n. mi) . A vertical angle of coverage limit Q was measured from 
the back-azimuth antenna. 

A new orthogonal coordinate system, with its origin located at the back- 
azimuth antenna and its X'-axis parallel to the runway center line, was used 
to make the geometric check. The north and east components of the estimated 

position vector Zg ^.-1 were transformed into the new X',Y’ coordinate system 
as follows: 

= Ze,N,t-l cos ^\)J. + Zg^E^t-1 sin 
Ye,t = Zg^N,t-l sin - Ze,E,t-1 cos 


The lateral azimuth check was then 


-Xg,t tan tJ ^ Yg,t = ^e,t tan w 
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The radial check was 


M ^ 


^e, t 
cos Tit 


^ L 


The vertical check was 


^a/p,t = *^o ^e,t ^ 


For these tests, the MLS back-azimuth coverage limits were 


Lateral (jJ = ±42° 

Vertical ^ = 20° 

Inner radial M = 3048 m 

Outer radial L = 37 040 m (20 n. mi.) 


RESULTS AND DISCUSSION 

Flight data were gathered on a total of 43 test runs along path 2 and 
path 3 and on the return flight to the start of path 1 . Position-estimate 
updates were accomplished on 27 update segments in the IBD mode, 7 update seg- 
ments in the IBX mode, and 30 update segments in the IDD mode. Typically, seg- 
ments lasted between 1 and 4 minutes. Position estimates frcan the navigation 
system in each update mode were compared with radar data to obtain the position- 
estimate error vector. The mean and standard deviation of the magnitude and 
direction of this vector were determined for each update mode. 


IBD Update Mode 

In the IBD mode, the magnitude of the position-estimate error vector had 
a mean of 124.3 m and a standard deviation of 60 m. The mean direction of the 
error vector was 298.1° with a standard deviation of 35.7°. 

Figure 8 shows how the position-estimate error vector is divided into two 
components, so that the errors due to the back-azimuth signal input and by the 
DME signal input can be separated. The first component ep is perpendicular 
to the back-azimuth radial on which the airplane is located. The ep component 
is influenced by the back-azimuth signal input. The second component ej^ is 
a radial error along the back-azimuth radial on which the airplane is located. 

The eR component is influenced primarily by the DME signal input and to a lesser 
extent by the back-azimuth signal. The eR component increases linearly as 
a function of the range R froii the back-azimuth antenna to the airplane if 
the azimuth bearing error e 0 (fig. 8 ) is constant (i.e., ep = Re 0 ) . 
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Figure 9 shows the ep component as a function of R for selected points 
from flights along the zero back-azimuth radial. These data are grouped about 
a straight line with a slope representing a constant e 0 = 0.05°. To determine 
how the ep component varied on other back-azimuth radials, data were used from 
flights along two other path segments. The first path segment was the last leg 
of path 3 which diverged frcan the runway center line at a 45° angle. The second 
path segment was parallel to the runway center line, but offset by 6100 m. This 
path segment was flown back toward the runway. These paths provided data at 
various ranges for all back-azimuth radials on the left-hand side of the azimuth 
coverage. (Flight during the test period was confined to this side of MLS cov- 
erage due to air-traffic-control considerations.) 

Figure 10 shows the ep component as a function of range from the back- 
azimuth antenna during flight along these path segments. The data lie above 
the line representing a constant 60 = 0.05°. Therefore, e 0 was larger for 
azimuth radials which were not on the runway center line. 

To determine how position error varied for different back-azimuth radials, 
the ep component from all flights (including those along the zero-azimuth 
radial) was normalized for range effects and plotted in the form of gq as a 
function of the back-azimuth radial. Figure 11 shows that these data were 
grouped about a straight line drawn from e 0 = 0.05° at the zero-azimuth 
radial to e 0 = 0.28° at the 40° radial. The straight line in figure 11 shows 
the expected e 0 error as a function of the back-azimuth radial. An expected 
ep component can be found by multiplying the range from the back-azimuth 
antenna and the appropriate e 0 obtained from the straight line. 

The radial component of error e^ was plotted as a function of range 
from the back-azimuth antenna in figure 12. The e^ components for flights 
on both path 2 and path 3 are randomly grouped about a mean value of 120 m. 

The figure shows that the e^ component did not vary significantly with range. 

The ej^ component contributed the most to the navigation position- 
estimate error. The ep component contributed a smaller error that varied 
as a function of the back-azimuth radial and the range from the back-azimuth 
antenna. The expected ep component was shown to be less than 32 m within 
a range of 37 040 m (20 n. mi.) from the back-azimuth antenna. For operations 
off the zero-azimuth radial, the ep component will become as large as the 
e^ component only vrtien the airplane is close to the lateral limits of the 
back-azimuth coverage and at large distances from the back-azimuth antenna. 

The DME ground station used during these flight tests was a standard 
EME station used for normal en route navigation. The expected accuracy of 
these stations is about 185 m (ref. 6 ) . A reduction of the EME ground- 
station error, which could be obtained with a PDME, could reduce the radial 
component of error and result in an improved navigation position estimate dur- 
ing the IBD update mode. 
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IBX Update Mode 

In the IBX update mode, the magnitude of the position-estimate error 
vector had a mean of 145.3 m and a standard deviation of 71.3 m. The mean 
and standard deviation for the direction of the position-estimate error 
vector was not summarized since radio position updating does not occur 
omnidirectionally in this mode. 

Both the magnitude and direction of the position-estimate error vector 
were strongly dependent upon the initial error when the IBX update mode 
started. Corrections to the position estimate were made only in a direction 
perpendicular to the back-azimuth radial on which the airplane was located. 

No corrections could be made in a radial direction since EME signal inputs 
were not available. 

Figure 13 shows the ep and ej^ error components for a flight along the 
zero-azimuth radial. The IBX update mode was started with initial navigation 
position-estimate errors of e^ = 140 m and ej^ = 118 m. The ep component 
was reduced because of back-azimuth signal inputs but the e^ component 
remained virtually unchanged since no DME signal inputs were available. 

Figure 14 shows the ep component as a function of range when the air- 
plane was flown along path 5 and back to the start of path 1 . The results 
were similar to those obtained in the IBD update mode: e 0 = 0.05° on the 
zero-azimuth bearing and e 0 ^ 0.25° on the return leg. 


IDD Update Mode 

In the IDD update mode, the position-estimate error vector was random 
in direction and had a mean magnitude of 183.5 m and a standard deviation of 
81.1 m. These results are within the 1 85-m accuracy tolerance for DME ground 
stations and are similar to past experience with the dual DME update mode 
(ref. 5) . It is expected that the use of PDME inputs also could improve the 
position-estimate accuracy in the IDD update mode. 


CONCLUDING REMARKS 

Flight tests with the NASA Terminal Configured Vehicle (TCV) Boeing 
737-100 airplane in the time-referenced scanning-beam (TRSB) microwave 
landing-system (MLS) environment have shown that the accuracy of position esti- 
mates based on dual distance-measuring-equipment (DME) inputs can be improved 
by using MLS back-azimuth radio signals. Position estimates from a combination 
of back-azimuth and DME signals achieved the best position-estimate accuracy. 
The next best accuracy was obtained from back-azimuth signals alone. The least 
accurate update mode of the three tested was the dual EME mode. 

The position-estimate accuracy was established by computing an error 
vector between the position of the airplane defined by ground-based radar 
tracking and the airplane position estimated by the airplane navigation system. 
This error vector was divided into two components: a component perpendicular 
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to the back-azimuth bearing, influenced by back-azimuth signal inputs; and a 
component parallel to the back-azimuth bearings, influenced primarily by the 
DME signal inputs. 

For updates using both DME and back-azimuth signal inputs, the radial 
component of error contributed the largest portion to the error vector. The 
radial ccmiponent of error was generally constant as a function of range from 
the back-azimuth antenna and was within accuracy tolerances for this DME 
station. The DME ground station used during these tests was a standard 
en route navigation station. The use of a precision DME ground station 
probably could improve the accuracy of the position estimate. 

The perpendicular component of error was influenced both by range from 
the back-azimuth antenna and by the particular back-azimuth radial on which 
the airplane was located. 

For updates which used only the back-azimuth signal inputs, position- 
estimate error was strongly dependent upon the initial position- estimate 
error since corrections were made only in the perpendicular direction. The 
perpendicular component of error showed accuracy characteristics similar to 
those observed during updates using both DME and back-azimuth inputs. 

For updates which used dual DME, position-estimate error was within the 
accuracy tolerance specified for en route DME and was similar to previous 
experience with dual DME updating. The accuracy of the position estimate 
could be improved with the use of precision DME signal inputs. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
Novonber 13, 1979 
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Figure 2.- Functional block diagram of navigation, guidance, and control systems. 
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Figure 3.- Flight test paths within the MLS front- and back-azimuth signal coverage. 

Headings are relative to true north. 
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Figure 4.- Geometry pertaining to IBD update mode: back-azimuth antenna, DME, estimated and measured 

positions of airplane, and runway with extended center line. 


















Figure 10.- Perpendicular component of position estimate error on various back-azimuth 

radials up to 40°; IBD update mode. 
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Figure 13.- The ep and e^ components of position-estimate error from an 

initial error; IBX update mode. 
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